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I. INTRODUCTION 


The general objective of this project is to advance our theoretical understanding of gas tori for 
the outer planet satellites and hence to enhance our ability to interpret observational data and to make 
available key concepts or quantitative results that are relevant to a number of other related theoretical 
studies. These studies are important to connect the aeronomy and photo- and ion chemistry of satellite 
atmospheres and their predicted gas loss rates, the evolution and distributions of gases in the larger 
planetary environment, and the ion loading and other impacts of these gases on the composition, 
structure, and properties of the planetary magnetospheres. Central objectives of this project are to 
explore for the satellites Titan and Triton the effects of two new mechanisms that we have very recently 
discovered (Smyth and Marconi 1993) to operate on gas tori in the outer planet satellite systems 
These two mechanisms can dramatically alter the current picture that has been widely adopted for the 
structure and evolution of long-lived gas tori (i.e., lifetime long enough to achieve approximate 
azimuthal symmetry about the planet) and. in particular, have profound consequences on the 
interpretation of Voyager data for both neutral and ionized species in the circumplanetary environments 
of Saturn and Neptune. Due, however, to the substantially reduced budget available for this project, 
the primary emphasis of the research was focused on the Titan component, with a limited effort (as 
time permitted) expended on the Triton component. 


For the Saturn system, comparable amounts of thermal H and H 2 and a much smaller amount 
of nonthermal N are thought to escape Titan and form gas tori. These tori are thought to be essentially 
collisionless because the lifetime loss processes of these gases in the moderate Saturn magnetosphere, 
although slow, are sufficiently rapid to avert collisional conditions. For the atomic hydrogen torus the 
perturbation of solar radiation pressure experienced by H atoms as they resonance scatter Lyma’n-oc 
photons - the first new mechanism - was shown by Smyth and Marconi (1993) to be operative and 
to destroy the normally assumed cylindrical symmetry of the torus produced by the 1/r central potential 
of a planet. This new mechanism causes H atom orbits to evolve inward as their eccentricities 
increase, and a significant fraction of these atoms is lost from the torus (having a preferred orientation 
of their perigee-axes) by colliding with the planet near its dusk side before they are otherwise lost 
through lifetime processes. This time evolution for a typical H atom lost from Titan is shown in 
Figure 1. Solar radiation pressure thus provides a natural mechanism for understanding the 
asymmetric distribution of hydrogen about Saturn recently reported by Shemansky and Hall (1992). 
This new understanding will subsequently allow the current question of the consistency of this atomic 
hydrogen distribution with the composition and properties of the inner magnetospheric plasma and the 
inner icy satellite ton to be addressed more clearly. In contrast to atomic hydrogen, H 2 and N will 

escape Titan and form approximately azimuthally symmetric tori about Saturn (i.e., the traditionally 


1 



adopted torus picture) since, for both species, solar radiation pressure is not important and the loss by 
lifetime processes is sufficiently long to achieve symmetry. The study in this project has been limited 

to the spatial nature the H gas tori of Titan to explore the consequences of the new "solar radiation 
pressure mechanism" noted above. 


For the Neptune system, comparable amounts of H and H 2 and N are thought to escape 
thermally from Triton and form gas tori. The lifetimes of these three gases in the magnetosphere of 
Neptune are, however, very long and are comparable to their photo-lifetimes because of the low 
plasma density and the small fraction of time that the tori spend in the more dense regions of the 
magnetosphere which executes complex motion about the planet. Because of these long lifetimes, the 
density of each torus will be colltstonal, and the neutral-neutral collision time will be shorter than all 
other time scales except for the typical Kepler orbit period of the atom or molecule (see Table 1 in 
Smyth and Marconi 1993). The collisional nature of the three gases will cause the multi-species gas 
torus to dynamically evolve in an inherently nonlinear manner that will depend upon the l/r nature of 
the potential of the planet. This collisional evolution will cause the gas torus to expand both inwardly 
(leading to inward gas loss by collision with the planet's atmosphere) and outwardly (leading to 
volume dilution and outward gas loss by escape from the planet). This expansion mechanism the 
second new mechanism noted above, is a non-linear effect that has been known for two decades in the 
field of solar system and ting formation but has been previously overlooked as important until now for 
the field of collisional gas tori about planets (Smyth and Marconi 1993). The final structure and 
density of the gas torus of Triton are particularly important in understanding the plasma sources 
recently identified in the analysis of the Voyager PLS data for the Neptune magnetosphere. This 
fascinating new expansion mechanism for the multi-species collisional gas torus of Triton provides a 
challenging dynamical evolution problem. Due to the reduced budget support available for this project, 
the solution of this problem to be undertaken in this project for Triton will be diminished in scope. 

e Triton component has been initiated only in this project with the limited resources focused on the 
initial phases of the development of a Direct Simulation Monte Carlo (DSMC) model to describe the 
collisional physics of the gas species in the circumplanetaiy volume. 

II. STUDIES FOR THE TITAN TORUS 


2.1 Model Description 

The basic model for the spatial distribution and Lyman-ct brightness of hydrogen in the Saturn 
system has been developed in the second and third project year. The model supports a central 
planetary source, an interior (combined inner icy satellite and ring) source, and a Titan source For all 
sources, a choice of a monoenergetic. Maxwell Boltzmann, or a sputtering initial velocity distribution 
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is possible. For the Saturn source, a sunlit hemispherical source is assumed and includes the spin rate 
and obliquity of the planet. For the composite interior source, a point source moving on a circular 
orbit about Saturn is used where the radial position of the orbit may be varied between 2 and 15 Saturn 
radii (R s ). For the Titan source, a point source moving on a circular orbit about Saturn of radius 
20.26 Rs is adopted, and the effect of the satellite's gravity on the initial exobase source velocity 
distribution is also included. The model calculates the density and brightness of hydrogen by 
following the time evolving trajectories of many H atom orbits initially emitted from the sources. The 
trajectories are calculated using the computationally efficient orbital element method discussed by 
Smyth and Marconi (1993) and dynamically include Saturn's spherical gravity, the leading 
nonspherical J 2 planetary component, and solar radiation pressure acceleration. Loss processes for the 
H atoms include reactions with the planetary magnetosphere and solar wind plasmas, collisional loss 
with Saturn and its rings, and escape of sufficiently energetic H atoms from the planetary system. The 
model of the Saturnian magnetosphere adopted in calculating the spatially nonuniform electron impact 

and charge exchange lifetimes for hydrogen is that of Richardson and Sittler (1990) and Richardson 
(1995). 

2.2 Observational Data for H 

The three best observations to date of the H distribution in the Saturn system were published 
by Shemansky and Hall (1992): (1) a Voyager 1 preencounter UVS scan providing an approximately 
in the satellite plane dawn to dusk Lyman-a brightness profile acquired in 1980 during the period day 
of year (DOY) 239 to 256, (2) a Voyager 2 preencounter UVS scan providing an approximately in the 
satellite plane dawn to dusk Lyman-a brightness profile acquired in 1981 during the period DOY 180 
to 186, and (3) a Voyager 1 postencounter UVS Lyman-a brightness mosaic image with the line of 
sight intersecting the satellite orbit plane at an angle of about 26“ looking down on the north pole of the 
planet acquired in 1980 during the period DOY 324 to 343. These observations show, in contrast to 
earlier reported Voyager data (Broadfoot et al. 1981; Sandel et al. 1982), that (1) there is no evidence 
for an absence of hydrogen in the radial interval inside of ~8 R s , but rather there is a peak in the 
Lyman-a intensity near the planet, and (2) there is no significant peak of Lyman-a intensity centered 
near Titan orbit as would be expected (Smyth 1981) based upon an azimuthally symmetric (doughnut) 
shaped hydrogen torus, but rather there is a local time asymmetiy in the Lyman-a intensity distribution 
with the preencounter scans showing a preponderance of emission on the dawnside of the system 
whereas the postencounter mosaic image data show more emission on the duskside. A reconciliation 
offered by Shemansky and Hall (1992) for the two different behaviors of the local time asymmetry 
was that the H distribution must have a complex three-dimensional morphology. Since at that time the 
satellite source was thought to produce an azimuthally symmetric H torus [the local time effect 
introduced by solar radiation acceleration (Smyth and Marconi 1993) had not yet been published], the 
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explanation for the nonuniformity in the local time distribution could only be attributed to a 
nonspherical Saturn source by Shemansky and Hall (1992). Model calculations undertaken in this 
project and discussed below which now include solar radiation acceleration, however, naturally 
produce a local time asymmetry in the circumplanetary H distribution for satellite sources and can fit 
the local time structures of the observations without a Saturn source. 

2.3 Model Calculations 

Modeling studies for the most definitive description of the spatial distribution of atomic 
hydrogen in the Saturn system discussed above and published by Shemansky and Hall (1992) have 
been undertaken and successfully completed in this project. These modeling studies are presented in a 
paper (Marconi and Smyth 1997) included in the Appendix. Only an overview summary of the major 
results and conclusions will therefore be presented below. 

The contribution of a Titan source to the preencounter UVS dawn to dusk scan Lyman-a 
brightness profile acquired by the Voyager 2 spacecraft in 1981 during the period DOY 180 to 186 is 
shown in Figure 1 by the solid line. An isotropic satellite exobase source strength of 4.8 x 10 27 H 
atoms s _l was assumed for a Maxwellian velocity distribution with an exobase temperature of 186 K. 
Because of solar radiation acceleration, the double-humped Titan source profile is asymmetric and 
higher on the dawn side (left side). The profile is able to fit the observed Lyman-a brightness at larger 
distances (beyond about 20 R s ), but is not able to fit the much brighter emission at smaller distances 
from the planet. An additional source of H atoms closer to Saturn and within Titan's orbit at 20.26 Rs 
is clearly required. The Lyman-a brightness profile, produced by combining an interior point source 
located on a circular orbit at 8.7 R s which corresponds to the orbital position of the icy inner satellite 
Rhea and the same Titan source adopted in Figure 1, is shown in Figure 2 and compared to the same 
Voyager 2 observation shown in Figure 1. The interior source is characterized by a 40,000 K 
Maxwellian velocity distribution, and its source strength is 1.9 x 10 28 H atoms s" 1 , a factor of about 4 
times larger than the Titan source. The contribution of the interior source is asymmetric and larger on 
the duskside of Saturn (right), opposite in since to the Titan source, and the combined interior and 
Titan sources fit the Voyager 2 observation in Figure 2 very well given the scatter in the data. 

A very good fit to the Voyager 1 preencounter UVS scan providing an approximately in the 
satellite plane dawn to dusk Lyman-a brightness profile acquired in 1980 during the period day of year 

(DOY) 239 to 256 is also achieved by using a combined interior source and a Titan source as shown in 
Figure 3. The best fit (solid line) was achieved for a slightly smaller Titan source of 3.3 x 10 27 H 
atoms s- 1 with a 186 K Maxwellian velocity distribution at the satellite exobase and a slightly smaller 
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interior source located on a circular orbit at 8.7 Rs with a source strength of 1.4 x 10 28 H atoms s _1 
still a factor of about 4 times larger than the Titan source. 

Assuming the same Titan source and interior source as in Figure 3, a model calculation for the 
Voyager 1 postencounter UVS Lyman-a brightness mosaic image with the line of sight intersecting the 
satellite orbit plane at an angle of about 26° looking down on the north pole of the planet acquired in 
1980 during the period DOY 324 to 343 has also been undertaken and is shown in Figure 4. Although 
there is substantial noise in the mosaic image, the only definite patterns evident is that both the model 
and observation are moderately brighter in the upper half of the ellipse defined by Titan's orbit (dotted 
line). It also appears that the overall magnitude of brightness variation is similar and that in both cases 
there appears to be a bright spot near the left hand apex of the satellite orbit. Within the limitations of 
the observations and the simplicity of the interior source adopted in the model, it would appear that the 
model and observation are broadly consistent. 

In order to test the effect of an interior source, a number of model runs were constructed with 
the interior source placed at different locations ranging from 3 to 15 R s . It was found that only if the 
source was placed in the neighborhood of 8 R s was a satisfactory fit to the observations in Figures 2 
and 3 achieved. This should not be interpreted as implying that the major source of hydrogen is in the 
neighborhood of Rhea’s orbital location. The precise distribution of the interior source is not well 
determined by the observations other than it must be largely confined inside of 8 or 9 R s . As for the 
interior source velocity distribution for hydrogen, a large number of different initial velocity 
distributions were tried including other Maxwellian with different temperatures as well as sputtering 
distributions. It was found that any distribution with a characteristic energy that was not too small (not 
much less than about 1 eV), could be combined with the canonical Titan source (used above) to 
produce a profile which was in reasonable agreement with the observations in Figures 2 and 3. The 

actual source of H is likely to have a complex velocity distribution since there are several different 
possible sources for hydrogen. 

The hydrogen source rates for the interior source above imply an H 2 0 production rate of -7-10 
xlO 22 molecules s-1 assuming H is created from H 2 0 and OH, with OH also produced mainly from 
H 2 0. If the dominant source of H is sputtering from icy surfaces with comparable amounts of H and 
H 2 also liberated, then this H 2 0 production rate could be reduced by a factor of 2 or 3 to a few xlO 27 
molecules s'l This lower H 2 0 production rate is still a factor of two or three higher than the updated 
maximum total sputtered rate given by Shi et al. (1995) for an interior source, although other sources 
such as micrometeoroids, not included in these estimates, may also contribute. The larger H 2 0 
production rate of -7-10 xlO 27 molecules s" 1 could be consistent with estimates for the atomic oxygen 


5 




loss rate of 3 xI027 atoms 5 -l within 4.5 R s of the planet by Shemansky and Hall (1992) if a 

comparable rate of loss occurs between 4.5 R s and the outer portion of our interior source at about 9 

Rs- The lower H 2 0 production rate of a few xlO^ molecules s‘l is. however, consistent with the 

amount of OH detected by Shemansky etal (1 993) at 4.5 R s as well as the amount of OH reported 
near Saturn's rings by Hall et al. (1996). 

Model calculations for a Saturn source of H were also undertaken with the result that Saturn is 
only a minor contributor at least outside of the immediate location of the planet where the error bars in 
the data m Figures 2 and 3 are rather large. The different local lime behaviots of the Triton source and 
the interior source produced by the acceleration of hydrogen by solar radiation pressure appear 
sufficient to reproduce the observations presently available. There is, however, current work in 

t ° fUrther redUCe lhese Voya * er observations for hydrogen in the Saturn system (Herbert 
1997 , and It is anticipated that this effort will eventually permit a more incisive comparison of the 
models and observations. 


II. STUDIES FOR THE TRITON TORUS 


The development of a model for the collisional evolution of gas tori in the Neptune system for 
given Triton sources was originally proposed to be undertaken by generalizing the numerical 
simulahou method of Trulsen (1972). In this project, however, an improved numerical approach has 

C T ‘ml T imPrOVed aPPr ° aCh brfcf1y diSCUSSed bel ° W “ Ca " ed DireCt Simulation Monte 
Carlo (DSMC as developed by Bird (1963), and effectively solves the non-linear collisional kinetic 

equation problem for a gas. Because of the limited resources available for this project, efforts have 
on T been able «o be focused on the initial phases of the development of a two-dimensional, multi- 
species DSMC model to describe the collisional physics of die gas species in the circumplanetary 
volume of Neptune. Progress in this development is summarized below. 

The determination of the evolution and steady state of the quasicollisional gas in the Titan torus 

DSMC h T Ve U " eXPl0red Pr0blOT ' ^ ,0 aCC ° mp,iSh ,hiS « have selected the 

DSMC approach which has a proven his.oiy of application to complex ramfied gas flow problems in 

he field of Aeronautics. In addition, contemporary DSMC techniques have attained a high degree of 

refinement and have been adapted to function on parallel computers. Because of these factors we 

ave opted to employ DSMC rather than the much less efficient and less general method of Trulsen 

We have been developing a spatially two-dimensional DSMC model for multi-species, chemically 

reacting, quasicollisional gases in a gravitational field. Since the Triton torus problem is highly 

demanding of computational resources, the model has been constructed to operate on parallel 
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computers so that the most powerful platforms can be brought to bear on the problem. At present the 
most basic elements of the model are completed and undergoing testing. In order, however, to apply 
this model to the Triton torus, several problem specific adaptations are necessary. These include the 
insertion of the Triton torus chemistry into the model and an appropriate Tritonogenic neutral gas 
source. Both of these tasks, although relatively straightforward, could not be completed with the 
limited resources available in this project and will be undertaken in a continuation of the Triton 
component of this project to be funded by the NSF Planetary Astronomy Program. 
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ABSTRACT 

The existence of atomic hydrogen in the Saturnian system has been known for 
almost a quarter century. Despite extensive observation and modeling over this 
period of time, the spatial morphology and kinetic nature of the circumplanetary 
hydrogen are still not well understood. The reasons for this lack of understanding 
are that the basic physical nature of the sources of hydrogen are poorly 
characterized and that the only detailed spatial measurements of hydrogen are 
from the Voyager excursions through the system. In this study, we combine the 
numerically efficient approach discussed in Smyth and Marconi (1993) for the 
calculation of orbits with the Saturnian plasma model of Richardson (1995) to 
model the Voyager 1 and 2 UVS observations as presented by Shemansky and 
Hall (1992). Although the comparison is limited by the quality of the data, source 
rates for a Titanogenic source of 3.3 - 4.8x1027 H atoms-s'l and, for the first time, 
source rates larger by about a factor of four for an interior source of hydrogen of 
1.4 - 1.9x1028 H atoms-s-l are obtained. These source rates for Titanogenic 
hydrogen are consistent with previous estimates. The source rates for the interior 
sources are consistent with recent estimates of source rates for other water group 
species based on observational consideration but are substantially higher than 
source strengths based on sputtering rate estimates. In addition to these source 
rates, some broad constraints on the energy and spatial distribution of hydrogen in 
the Saturnian system are also given. 



1. INTRODUCTION 


The existence of atomic hydrogen in the Saturn system was first predicted by 
McDonough and Brice (1973). Two years later, in 1975, rocket borne 
observations of Lyman-a emissions from near Saturn’s rings (Weiser et al. 1977) 
detected atomic hydrogen for the first time. In 1976 and 1977, the earth-orbiting 
Copernicus satellite (Barker et al. 1980) provided evidence for hydrogen near 
Titan, in 1979 the Pioneer 1 1 spacecraft found hydrogen near Titan's location on 
its orbit and near the rings (Judge et al. 1980), and from 1978 to 1990 (McGrath 
and Clarke 1992), the International Ultraviolet Explorer (IUE) has been utilized at 
least once per year to observe the hydrogen near Saturn’s disk. The most spatially 
detailed observations of hydrogen in the Saturnian system, however, remain those 
of the Voyager 1 and 2 missions that flew through the Saturn system in 1980 and 
1981. These observations provided one-dimensional scans of the Lyman-a 
emissions from hydrogen throughout the entire Saturnian system (Broadfoot et al. 
1981, Sandel et al. 1982) and also provided a detailed picture of the 
magnetospheric plasma (Bridge et al. 1981, 1982; Richardson and Sittler 1990) 
which is a critical component in understanding the sources and sinks for this 
hydrogen cloud distribution. As a result, most studies to date have focused upon 
the Voyager data in their attempts to explain the source and structure of this 
hydrogen. 

In conjunction with the foregoing observational efforts, there has been 
substantial attention directed to a theoretical understanding of hydrogen in the 
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Saturnian System. The earliest studies (McDonough and Brice 1973 and Fang et 
al. 1976, for example) assumed that the only source of hydrogen was Titan, that 
■he hydrogen was collisionless, that the only force acting on the hydrogen atoms 
was Saturn's gravity, and that the loss rate of hydrogen was, moreover, 
characterized by a single spatially homogeneous lifetime. Studies of Dennefeld 
(1974) considered a planetary source, a ring source, and a Titan source and 
concluded that the Saturn source was not significant, that an H 2 0 atmosphere 
would exist for the rings, and that an H and H 2 torus would exist for Titan. For a 
Titan source, Smyth (1981) also included the effect of the gravity of Titan on the 
hydrogen atom orbits and calculated in three-dimensions the time evolution of the 
torus, showing that the circumplanetary hydrogen distribution for an appropriate 
H lifetime was essentially azimuthally symmetric with a density maximum near 
T.tan's orbit. Generally speaking, these models led to a torus of hydrogen about 
Titan's orbit which was approximately axisymmetric (doughnut shaped), devoid 
of significant density mside of ~8 planetary radii, and whose density was 
determined by the source rate and assumed lifetime. This picture was broadly 
confirmed by the preliminary analysis ofBroadfoot et al. (1981, 1982) of some of 
the Voyager 1 and 2 ultraviolet spectrometer (UVS) observations. Their analysis 
concluded that the hydrogen in the Saturnian system formed a roughly unifotm 
torus, extending from 8 to 25 R s (Saturn radii), with a cavity within 8 R s , and a 
typ cal density of 20 cm 3. Although a post- Voyager model of the Titan H torus 
by Ip (1985) incorporated a one-dimensional (radial) lifetime profile, this picture 
of the hydrogen cloud dominated the general view until Shemansky and Hall 
(1992) analysed previously unpublished UVS observations and clearly 
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demonstrated that the hydrogen extended to Saturn and was characterized by a 
complex three dimensional structure quite unlike a doughnut. With this new 
asymmetric hydrogen cloud in mind, Smyth and Marconi (1993) studied the 
particle orbits of H atoms emanating from Titan. They introduced two new 
ingredients, namely, radiation pressure and the leading non-spherical component, 
J 2 , to Saturn's gravity. While they did not specifically calculate the hydrogen 
cloud, they noted that the orbits obtained in this scenario would lead to a highly 
asymmetric cloud extending to Saturn. Ip (1995) explicitly calculated the 
hydrogen cloud produced by a Titan source for the scenario in Smyth and 
Marconi and confirmed that the cloud was quite asymmetric and extended as far in 
as Saturn. Ip also compared the radial density profile of the cloud to the UVS 
data presented by Shemansky and Hall and concluded that another source of 
hydrogen was necessary within 1 5 R s . 

Over the same period, other studies suggested that Titan was not the only 
source of atomic hydrogen in the Saturnian System. Shemansky et al. (1985) 
analysed Voyager 2 data and showed that atomic hydrogen extended into Saturn's 
atmosphere and Hilton and Hunten (1988) applied their model for hydrogen to the 
Voyager measurements and concluded that in order to explain the vertical extent of 
the hydrogen, Saturn must be the principal source of atomic hydrogen within 10 
R s . Further indirect support and a mechanism for Saturn as a major source of 
atomic hydrogen came from the Lyman-oc dayglow phenomena observed at 
Saturn. Failure to account for the dayside Lyman-a dayglow by scattering of 
solar Lyman-a radiation led to the suggestion that collisional excitation of atomic 
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hydrogen by electrons in the upper atmosphere of Saturn (electroglow) could 
account for the excess Lyman-a brightness of Saturn (Shemanksy and Ajello 
1983, Shemansky 1985, Shemansky et al. 1985, and Yelle et al. 1986). 
Furthermore, these electrons could produce, by impact dissociating H 2 , enough 
sufficiently energetic hydrogen atoms to account for the atomic hydrogen cloud 
within -8 R s - Recently, however, Emerich et al. (1993) applied a new and 
improved radiative transfer model for the scattering of solar Lyman-a by H to 
various Lyman-a airglow data for the giant planets and concluded that the need 
for the electroglow process was much reduced. Ben Jaffel et al. (1995), using the 
same approach in a new analysis of the Voyager UVS H Lyman-a emission from 
Saturn, concluded that reflected Lyman-a was the main contributor to the 
brightness rather than collisional excitation. Hence, a mechanism for dispersing 
hydrogen atoms from Saturn to large distances is not presently clear. 

The Voyager 1 and 2 plasma observations (Bridge et al. 1981, 1982) also 
focused attention on the surfaces in the Saturnian system as additional sources of 
atomic hydrogen. The rings of Saturn (Ip 1984) and icy surfaces of the Saturnian 
satellites, as a result of bombardment by magnetospheric plasma, solar radiation, 
and micrometeoroids, could directly generate atomic hydrogen or other heavier 
neutral species (Cheng et al. 1982, and Lanzerotti et al. 1983). Furthermore, gas 
phase chemical interactions between the heavier species, such as H 2 0, with the 
plasma and solar radiation were also recognized as another source of atomic 
hydrogen. Chemical "one box" models of the magnetospheric ions and neutrals 
were constructed by Richardson et al. (1986) and Richardson and Eviatar (1987) 
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who concluded that the neutral density had to be quite small and supported the 
view that Titan was the major source of hydrogen in the Saturn magnetosphere. 
Johnson et al. (1989) calculated the spatial distribution of plasma sputtered H 2 O 
from the icy Saturnian satellites and found high density neutral tori associated 
with the satellites. Shemansky and Hall (1992) using a one box chemical model of 
the ions and neutrals in the Saturnian magnetosphere, which included some key 
chemistry neglected in the Richardson models as well as a more reasonable 
estimate of the ion diffusion rate, concluded that the neutral densities were much 
greater than the Richardson models predicted. The high density of neutrals in the 
Saturnian magnetosphere was recently confirmed by the Hubble space telescope's 
faint spectrograph observation of the OH radical near the orbit of Thetys 
(Shemansky et al. 1993). The large abundance of OH inferred from this 
observation created a problem in that the estimates of neutral production rates 
from icy satellite sputtering, as calculated by Johnson et al. 1989 for example, 
were well below the level required to sustain it. Shi et al. (1995) used new, 
substantially larger, sputtering rates derived from some recent laboratory 
measurements to determine the sputtering rates of neutrals from the icy satellites. 
While the rates were considerably higher, they still fell quite short of explaining 
the OH observation. Most recently, further observations by Hall et al. (1996) 
with the Hubble Space Telescope detected in addition large amounts of OH above 
the A and B rings of Saturn. 

To summarize, it seems that there is general agreement that Titan is a 
substantial source of hydrogen for the Saturn magnetosphere. This hydrogen 
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probably dom, nates beyond roughly 10 R, and is generated a, a rate roughly in the 
range of 2-5x1 0 2, atoms-s'. In the inner magnetosphere, however, other sources 
are dominant. The nature and strength of these sources are far less certain. 
Saturn's atmosphere would seem to have an abundance of hydrogen to supply ,o 
• h = inner magnetosphere. „ is no, clear, however, that the electroglow process is 
able to provide the requisite flux of hydrogen and, moreover, the recent reanalysis, 
mentioned above, of dayglow observations do no, requtte a significant electroglow 
contribution to the Saturnian Lyman-o brightness any case. The icy satellites 
mgs are sputtering sources of neutrals for the magnetosphere of Saturn In 
.his case, there is .Me problem with populating the inner magnetosphere with 
hydrogen atoms either directly or from the destruction of hydrogen bearing 
spectes such as H j0 . Up to the present, however, the sputtering source rates 

required to explain the neutral observations are substantially above Are sputtering 
rates based on laboratory data. 


he purpose of this paper is to address the nature and strength of the Saturn 

nng, icy satellite, and Titan sources of hydrogen in the Saturnian magnetosphere 

A model which includes a Titan source, a Saturn source, and a source in die icy 

satellite region and ring region between Saturn and Titan (i.e„ an interior source) is 

utilized to obtain source rates and some general characteristics of the interior 

sources by using Voyager observations of Ihe hydrogen distribution enveloping 
the Saturnian system. 
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2. MODEL 


The basic model that is used in this work essentially computes the density and 
brightness distribution of hydrogen emitted from a source under the assumption 
that the hydrogen is to a good approximation collisionless (Smyth and Marconi 
1993). Following Smyth and Marconi (1993), the model includes Saturn's gravity, 
including the leading nonspherical, J 2 component, solar radiation pressure, loss of 
hydrogen atoms to reactions with magnetospheric and solar wind plasma, and 
collisional loss of hydrogen to Saturn and its rings. There are three basic steps to 
performing the model calculations presented below: (1) choosing the initial 
conditions for H atoms to be emitted, (2) calculating the orbits corresponding to 
these initial conditions, and (3) calculating orbital information so that spatial 
hydrogen densities and brightnesses may be determined. 

2.1 INITIAL CONDITIONS 

The first step in the model is to choose the initial velocities, positions, and 
emission times for the collection of hydrogen atoms, typically a hundred 
thousand, emitted from one of the three sources considered here, i.e. Titan, Saturn, 
or some interior source placed between Titan and Saturn. This process basically 
involves choosing the speed randomly from some speed distribution, for example 
Maxwellian, and the direction according to some angular distribution for the 
emission of a hydrogen atom. While various speed distributions were considered, 
the angular distribution was always taken as isotropic. In the case of Titan, the 
initial velocity distribution was based on an exospheric temperature of 186 K and 



an exobase radius of 41 75 km (Table I). Assuming an isolropic Maxwellian with 
a temperature of 186 K at the exobase, the resulting veloeity distribution at the 
Lagrange radius was approximated by reducing the energy of the particle by an 
amount equal to the gravitational potential energy between the exobase and 
Lagrange radius. A separate three body (including the gravity of Saturn and Titan) 
more exact calculation for the orbits between the exosphere and Lagrange radius of 
Tttan, assuming a 186 K Maxwellian at the exobase of Titan, was also performed 
and confirmed that the simple reduction in energy yields a good approximation to 
the actual velocity distribution at the Lagrange radius for this case. The Titan 
source was then treated as a point in these calculations and the point of emission 
was the position of Titan in its orb,,. Emission times were chosen at random over 
a time interval chosen sufficiently long so that initial trans.ents would decay 

(~5xI0 s). This was well satisfied for the chosen emission time interval of 2x10* 
s for particles emanating from Titan. 

The interior source was treated almost identically to Titan. It was also taken 

in the planetary equator plane as a point source moving in a cireular orbit around 

Saturn. The radial position of the source was varied between 2 R, and 15 R,, 

however. In addition, the initial speed distributions of the emitted hydrogen, 

which included Maxwellians and modified sputtering distributions (Smyth and 

Combi 1988), were varied in contras, to the 186 K Maxwellian always assumed 
for Titan. 
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Finally, Saturn was treated as a source with radius, spin rate, and obliquity 
given in Table 1 . The Saturn source was assumed to be emitting preferentially 
from the sunward hemisphere as suggested by observations (Shemansky and Hall 
1992, for example). The source strength dependence on location on the sunward 
hemisphere was assumed uniform due to the lack of any position dependent 
information and also due to a lack of sensitivity of the results to the spatial 
dependence of the hydrogen source. The initial hydrogen speed distributions 
were, as for the interior source, varied. Since satellite orbits are difficult to 
populate (Hilton and Hunten 1988) only escape and ballistic orbits were 
considered. The effect of the rings was ignored since only a small fraction of 
orbits intersect the rings. Moreover, of the small fraction that strikes the rings, 
much of that will be scattered back towards Saturn which fills a large part of the 
sky as seen from the rings. 

2.2 ORBIT CALCULATION 

Once the initial velocity, position, and time of emission are determined, as 
described above, the orbit integration is performed. In this work, rather than 
following the spatial position of the atoms by directly time integrating the 
equations of motion, the initial conditions are converted into orbital elements and 
the evolution of the orbital elements in time are integrated as in Smyth and 
Marconi (1993). This approach accounts for the perturbative effects due to both 
the J 2 component of gravity and radiation pressure (values give in Table 1) which 
were shown to be important by Smyth and Marconi (1993). The advantage over 
the conventional approach of solving directly Newton's equations for the orbits is 
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that the integration proceeds much more quickly and, thus, permits a much larger 
number of particles to be used and many more scenarios to be tested. 

2.3 DETERMINATION OF DENSITIES AND BRIGHTNESSES 

In order to determine the hydrogen density or brightness, the loss of hydrogen 
must be calculated as the atoms travel along their orbits. As mentioned above, the 
model accounts for the loss of hydrogen from reactions with the magnetospheric 
and solar wind plasmas, collisions with Saturn and its rings, and escape of 
sufficiently energetic hydrogens from the system. 

The model of the Saturnian magnetosphere used in this study is that of 
Richardson and Staler (1990) and Richardson (1995). The equatorial densities, 
temperatures, and velocit.es, of the heavy ions, protons, hot electrons, and cold 
ctrons from L shell, L - I to L = 20, were taken from Richardson (1995). 
These quantities were then used with Richardson and Sittler-s (1990) model of 
force balance along the magnetic field line to genetate the plasma properties above 
and below the equator plane. This field of plasma ptoperties was converted into a 
field of loss rates (see Figure 1) which include charge exchange with the heavy 
tons, assumed to be O’, charge exchange with protons, as well as impact 
ionizatton by the hot and cold electrons. Richardson (1995) provides two 
alternative models for the equatorial plasma properties between L = 12 and L = 

20. Since model calculations turned out to be similar in both cases, only one 

(Richardson's alternative 1) of these alternative models is used in the 
computations below and in Figure 1. 
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The plasma information used here is derived from the Voyager plasma science 
instrument (PLS) and is restricted to an energy/charge range of 10 to 5950 Volts 
(Richardson and Sittler, 1990). As a result, the true loss rates should be higher 
since there is plasma with energy outside of the PLS energy window. The amount 
of plasma flux outside of the PLS window should be small, however, judging by 
the flat flux energy spectra toward the low energy limit of PLS, the rapid drop 
toward the high energy limit of PLS, and the much lower value of plasma fluxes 
measured by the Voyager low energy charged particle instrument (LECP) starting 
at -100 keV (Shi et al. 1995). As a result, the magnetospheric plasma loss rate of 
hydrogen from plasma within the PLS energy window is probably a good 
approximation to the total loss rate unless there are unexpectedly large plasma 
fluxes at the unviewed energies. 

Outside of the spatial range of the above PLS data, i.e. magnetic flux shell L = 
20, the loss rate of hydrogen was assumed to be due to resonant charge exchange 
with solar wind protons and a nominal value of 2x10' s' was most often 
employed. This choice for the loss rate is similar to theoretical estimates of the 
rate assuming nominal solar wind parameters at earth of 6 protons cm' 3 , a velocity 
of 400 km-s' 1 , a charge exchange cross-section (Stebbings et al. 1964) of 7.5xl0' ls 
cm' 2 , and then scaling to Saturn's heliocentric distance. This is also similar to the 
rate used by Ip (1995) of 3x10' s' beyond 25 R*. Some models were also 
constructed for different values of the loss rate. The best agreement with data 
was, however, obtained for a loss rate near 2x10’ s' . 
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Because, for example, of the changing location of the solar wind bow shock, 
the plasma beyond 20 R s (Bridge at al. 1981, 1982) is quite variable in space and 
time so that it is difficult to accurately construct a suitable representation of the 
plasma which could be utilized in this model. It appears, however, that a broad 
spatial and temporal average of the measured plasma properties leads to a loss rate 
of the order of 10' 8 s '. Moreover, the large spatial extent and long lifetime of the 
orbits penetrating this region means that only the average properties of the plasma 
are important rather than the smaller scale spatial and faster temporal plasma 
variations. Consequently, the choice of a spatially and temporally uniform loss 
rate with the above value is well justified. For all orbits, the relatively small 
contribution of solar photoionization of the hydrogen with a rate of 7.26x1 O' 10 s -l 
(Huebner et al. 1992) was also included. 


The actual determination of the amount of hydrogen lost from reactions with 
the plasma is performed in three steps. In the first step, the orbital elements at 
the start of some integration time step are used to calculate the spatial location of 
the corresponding Kepler orbit. The orbit is then divided into a large number of 
uniform time intervals and the location of each time interval in space is computed. 
The small loss of hydrogen in each of these time intervals is determined with the 
loss rate model discussed above and then summed over each time interval to g,ve a 
total loss for the orbit. In the second step, the total loss over an integration time 
step is calculated. The orbital elements are assumed to be unchanged over an 
integration time step so that the same Kepler orbit may be used. The total loss 
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over this integration timestep is then obtained by suitably scaling the orbital 
period with the ratio of integration timestep to the period of the Kepler orbit. In 
the third step, this loss is added to the amount of hydrogen present at the end of 
the previous time step and the procedure restarted with new updated orbital 
elements appropriate to the time corresponding to the start of the next integration 
time step. 

Additional losses of hydrogen occur due to collisions with Saturn and its rings. 
Collisions with Saturn are instantly evident, when in the course of integrating the 
orbital elements, the eccentricity of the orbit becomes sufficiently close to 1. In 
this case, the particle is assumed to be absorbed by Saturn. In the event that the 
orbit of a particle intersects the rings of Saturn, the procedure followed is similar 
to Ip (1995), i.e. if a particle crosses the equator within 2.27 R* the probability of 
absorption is 0.26 corresponding to an optical depth of 0.3. Since the results are 
not particularly sensitive to the presence of rings, this simple model of ring 
absorption is sufficient for this study. Finally, if the atom is on an escape orbit, 
the orbit is tracked to some distance slightly greater than the region of interest 
(±40 Rj from Saturn) and then the atom is assumed lost. 

In order to calculate density and Lyman-a brightness, a three-dimensional grid 
is overlaid on the collection of orbits generated in the orbit integration part of the 
model described above. The grid used in this study consists of 80 x 80 x 60 boxes 
in the x, y, and z directions of size 1 R^x 1 RjX 1 Rs respectively. The spatial . 
coordinates (x,y,z) were chosen to define an inertial coordinate system centered at 
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Saturn with x and y in the equatorial plane and z perpendicular. The hydrogen 
density in this grid is determined by recording whether an orbit traverses a 
particular box and then correcting for the fraction of hydrogen remaining at that 
point due to lifetime processing. The contribution per unit volume to the 
brightness, i.e. the amount of solar Lyman-a resonantly scattered, is determined 
by multiplying the density by a g-factor which accounts for the heliocentric 
distance of Saturn and relative velocity of the atom and the Sun. This is similar to 
Smyth et al. (1995) except that the g-factor used here was normalized to the time 
dependent solar Lyman-a flux averaged over the observation period and scaled to 
Saturn's heliocentric distance (Hall 1997 and Pryor 1997). The contributions per 
unit volume to the brightnesses in the grid boxes, described above, are summed 
along the appropriate line of sights and convolved with the appropriate Voyager 
slits, given in Shemansky and Hall (1992), to finally obtain the brightnesses for 
the images and radial profiles shown below. 

Before we proceed to a presentation of the results of this study, there are two 
items which should be mentioned. The first is that impinging on the Saturnian 
system there is a flow of interstellar hydrogen which could effect the distribution 
of hydrogen emanating from the sources of hydrogen in the Saturnian system. 
Using nominal values from Allen (1976) of a density of 0.1 cm* 3 and velocity of 
21 km-s for the interstellar hydrogen, the collision rate with the Saturnian 
system hydrogen is of the same order of size as the photoionization rate, i.e. 

10 s . This may be neglected in comparison to the smallest rate in our plasma 
loss rate model of 2x1 0‘ 8 s' 1 for hydrogen in the solar wind plasma. 
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Another potentially important issue regards the neglect of the effect of Titan 
on the hydrogen orbits. A hydrogen wandering sufficiently close to Titan could 
be captured or strongly scattered. If we assume that this "sweeping" radius of 
Titan is equal to its Lagrange radius, ~20 Titan radii (Rt) (Table I), then model 
computations of the amount of hydrogen in the neighborhood of Titan's orbital 
location indicate that if it is assumed that this hydrogen is stationary in the 
Saturnian inertial frame. Titan's Lagrange sphere would intercept an amount of 
hydrogen equivalent to roughly an amount comparable to the hydrogen emanating 
from Titan. This would have significant consequences for the hydrogen cloud. In 
reality, however, most hydrogen atoms have sufficient velocity relative to Titan 
so that they would have to penetrate to well within the Lagrange sphere before 
being substantially perturbed by Titan's gravity. If we assume that the mean 
speed, relative to Titan, of a hydrogen atom in the neighborhood of Titan's orbit is 
about the same as the mean speed of the initial distribution of hydrogen atoms 
from Titan at the Lagrange radius of ~1.4 km-s -1 , then a hydrogen atom would 
have to penetrate to ~4 to 6 Rj from Titan before being significantly disturbed by 
Titan's gravity. Therefore, it is likely that the effective sweeping rate of Titan is a 
small fraction of the production rate of hydrogen from Titan and the neglect of 
Titan's gravity in this work is justified. 



3. RESULTS 


Figure 2 is a comparison of the Lyman-a model brightness (solid line) of the 
hydrogen cloud due to a purely Titanogenic source of 4.8 x 10 27 atoms s' 1 with the 
observed brightness (crosses) obtained by accumulating the Voyager 2 UVS 
approximately radial scans taken from day of year (DOY) 180 to DOY 186 in 
1981 as described in Shemansky and Hall (1992). The brightness asymmetry in 
the double-peaked model profile is produced by the action of solar radiation 
acceleration. Clearly, the match is inadequate and strongly suggests that an 
additional, more centrally located source of hydrogen is required, as was also 
pointed out by Ip (1995). 


In Figure 3, the same data is the subject of comparison. The model, however, 
now includes two sources. The same Titanogenic source that was used in Figure 2 
is combined with an internal source of 1.9xl0 28 atoms s* 1 at 8.7 R s which 
corresponds to Rhea's orbital position. This interior source is assumed to generate 
hydrogen atoms with a 40000 K Maxwellian. The matching of the model to the 
data displayed in Figure 3 is achieved by adjusting the ratio of the Titan to interior 
source contribution and overall normalization until the model brightness profile 
approximates some mean trajectory through the region defined by the data points 
and their error bars. The data profile displays significant oscillations which are 
not evident in the model. These oscillations, however, are largely spurious 
(Doyle, 1996) and a mean curve through the data is an appropriate basis for a 
comparison. In Figure 4, the same two sources are similarly adjusted with slightly 
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smaller source rates to match the nearly radial UVS scans accumulated to produce 
the H Lyman-oc intensities obtained by Voyager I during the period DOY 239 to 
256 in 1980 (Shemansky and Hall, 1992). For both sets of data. Figure 3 and 
Figure 4, the agreement of the model with the data is quite good and certainly is 
vastly improved over the case of the Titanogenic source alone shown in Figure 2. 

The choice of the 40000 K Maxwellian for the interior hydrogen source 
represents the best of a large number of other models with different initial velocity 
distributions. A wide variety of other velocity distributions were tried which 
included other Maxwellians with different temperatures as well as various 
sputtering distributions. It was found that any distribution with a characteristic 
energy that was not too small, i.e. much below about 0.5 eV, could be combined 
with the canonical Titan source, described above, to produce a profile which was 
in reasonable agreement with the data shown in Figures 3 and 4. In addition, the 
source rates inferred for the other distributions were similar to those obtained for 
the 40000 K case. The actual source of H is likely to have a complex velocity 
distribution since there are various different sources for hydrogen. H may be 
generated from H 2 0, OH, and H 2 as a result of photodissociation or electron 
impact dissociation. H may also be directly generated from surface sputtering. 
All these sources contribute H atoms with velocity distributions characteristic of 
the particular source. Moreover, since H is light, it will likely not be spatially 
localized near its source. As a result, the determination of the nature and 
distribution of the sources of H is not a straightforward task, particularly in view 
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of the limited sensitivity of the data constdered here to the details of the velocity 
and source distribution for H. 


In order to test the effect of the placement of the interior source, a number of 
models were constructed with the interior source placed at different locations 
ranging from 3 R, to 15 R,. It was found that only if the source was in the 
neighborhood of 8 R,, as in the interior model used in Figures 3 and 4, was a 
satisfactory match obtained. This, however, should no, be interpreted as 
implying that the major source of hydrogen is in the neighborhood of Rhea's 
orbital location. As for the velocity distribution of hydrogen, there is limited 
sensitivity to the choice of source. The exact spatial distribution of the source is 
not well determined by the data model comparison displayed in Figures 3 and 4. 
The obtaining of a better match with an interior source at 8.7 R, actually only 
indicates that the source of non-Titanogenic hydrogen is hugely confined to inside 
of around 8 or 9 R,. The precise distribution of the source within that radial 
distance is not well determined in this comparison. Models with a Saturn source 
of H were also calculated with die result that Saturn is only a minor contributor a, 
leas, outside of die immediate location of Saturn where the etror bant in the data 
are rather large. This is consistent with the Lyman-a dayglow being due to 

scattering of Lyman-a radiation rather than collisions processes which could 
produce energetic hydrogen. 


Numerical models were also computed with different loss rates for H in the 
solar wind plasma. A definite preference for a loss rate 1.5 - 2.0x1 O’ 8 s' 1 (lifetime 
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5xl0 7 s) was noted. Other choices of loss rates, significantly outside this range, 
yielded clearly inferior agreement with the data in Figures 3 and 4. It is interesting 
to note that this rate is very similar to the simple scaling to Saturn's heliocentric 
distance of the charge exchange rate for H near the Earth using nominal solar wind 
conditions. 

Source rates for H obtained for Titan and the interior source were, 
respectively, for the data in Figure 3, 4.8xl0 27 atoms s’ 1 and I.9xI0 28 atoms s' 1 , 
while for the data in Figure 4 were, respectively 3.3xl0 27 atoms s’ 1 and 1.4xl0 28 
atoms s' 1 . In both cases the interior source is about 4 times larger than the Saturn 
source. This discrepancy (-40%) in the absolute source rates between the two 
dates could arise from a number of factors associated with the modeling such as 
the matching of the model to data which is an estimation by eye of a reasonable 
approximation to the highly oscillatory data. In addition, the slit geometry is 
complicated with only representative examples shown in Shemansky and Hall 
(1992). It also appears that the data processing performed by Shemansky and 
Hall is quite involved so that together with the modeling performed here, there is 
ample opportunity for spurious variations to enter the data model comparisons. 
In principle, since the measurements were well separated in time, intrinsic 
variation of the source could be also possibility. It is not clear why both sources, 
however, should be affected and in almost the same proportion. Another 
possibility is that the solar Lyman-a brightness information used in this study is 
based on the He 10830 proxy method rather than direct measurement of solar 
Lyman-a luminosity. Since the Lyman-a solar luminosity is known to vary by 
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as much as 40% over a solar cycle (Skinner et al. 1988), this represents another 
potential source for the differences in the inferred source rates for H. 

The source rates derived for the Titan source are in the range of 2-5x1 0 27 
atoms s 1 which, as mentioned above, are in the generally expected range. The 
hydrogen source rates for the interior source correspond to an H 2 0 source rate of 
7-10x10 molecules s 1 assuming H is mainly from H 2 0 and OH (assuming OH is 
also mainly from H 2 0). This is an order of magnitude higher than the maximum 
total sputtering rate given by Shi et al. (1995). According to Bar-Nun et al. 
(1985), however, comparable amounts of H, and H 2 0 may be sputtered as well. 
If this is true, the actual amount of H 2 0 required may be reduced by a factor of 2 
or 3 to a few x 10 molecules-s 1 so that the discrepancy with the estimates of Shi 
et al. may be closer to a factor of 2 or 3. While Shi et al. have employed the latest 
laboratory data to estimate these rates, we repeat the comment by Shemansky and 
Hall (1992) that there is an inherent difficulty in translating laboratory surface 
reactions to icy satellite surfaces. In addition, there are other means of ejecting 
H 2 0 molecules from the surfaces in this system, such as micrometeoroids, not 
treated in their work. The largest H 2 0 production rate of 7- 10x1 0 27 molecules s' 1 
estimated above is, however, consistent with Shemansky and Hall (1992) who 
estimated a loss rate of O of 3xl0 27 atoms-s' 1 within 4.5 R s , if a comparable rate of 
loss occurs between 4.5 R, and the outer position of our interior source at about 9 
Rs- It appears that larger rates of H 2 0 production than can be accounted by Shi et 
al. are also implicated in the amount of OH detected by Shemansky et al. (1993) at 
4.5 R s as well as the amount of OH reported near Saturn's rings by Hall et al. 
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(1996) as concluded by Ip (1997). In both cases, a factor of 2 or 3 more 
sputtering of H 2 0 is required which is consistent with our findings if significant 
amounts of H and H 2 are also sputtered. 

Figure 5 is a contour plot of the relative solar Lyman-a brightness 
corresponding to the observing geometry for the image displayed in Figure 1 of 
Shemansky and Hall, i.e. during the period 1980 DOY 324-343. The solid circle in 
the center is Saturn and the dotted line is Titan's orbit as would appear from 
Voyager's vantage point. The Sun is in the direction of the arrow. The same 
mixture of Titan source and interior source was used as in Figure 4. This 
comparison of the model with image is highly desirable as it would ordinarily 
represent a much more severe test of the validity of the model than the one 
dimensional data model comparisons discussed above. In this case, however, there 
appears to be substantial noise in the image. The only definite pattern evident is 
that both model and data are moderately brighter in the upper half of the ellipse 
defined by Titan's orbit (dotted line in Figure 5). It also appears that the overall 
magnitude of brightness variation is similar and that in both cases there appears to 
be a bright spot near the left hand apex of the orbit. Within the limitations of the 
data and the simplicity of the sources in the model it would appear that the model 
and data are broadly consistent. 

The enhanced brightness in the region above Saturn in Figure 5 also suggests 
that the interior source is much stronger than the Titanogenic source. The orbital 
dynamics of an atom produced by gravity and radiation pressure leads to an 
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orb.tal evolution as shown in Figure 6 (see also Smyth and Marconi, 1993) which 

is for a typical Titanogenic hydrogen. Figure 6 it is seen that there is an 

accumulation of material on the left of the Sun-plane, line. A large fraction of this 

mater, al is eventually lost by collisions w„h Saturn. Since the temperature of the 

Titanogenic source is relatively small, the hydrogen atoms from Titan tend to 

execute orbits with small eccentricities and thus avoid the high density plasma 

•nside of Titan's orbit. As a result, the plasma lifetime of the Titanogenic 

hydrogens is eomparable to the time requited for an orbit to undergo the evolution 

depicted Figure 6 which requires ~,0> s. The result is an accumulation of 

materia, to the left of the sun plane, line as implied in Figure 6. This Titanogenic 

hydrogen corresponds below Saturn on Figure 5 and is contraiy to what is 

observed. On the other hand, the orbits generated from the interior source are of 

high eccentricity due to the assumed high temperature of the interior relatively 

source and thus pass through the dense p,asma where the hydrogen lifetime is 

short (Figure 1). While these orbits engage in the same evolution as depicted in 

Ftgure 6, they are more quickly destroyed by the plasma so tha, a, any given time 

the distribution of surviving orbits have evolved for a time short compared to the 

olution time. As a result, there will be more interior source hydrogen atom 

orbits to the right of the sun planet line since some of the orbits to the left will 

have been lost by collisions with Saturn Thus the in^r 

• tnus, the interior source tends to 

produce a somewhat larger brightness above Saturn as seen in Figure 5 This 
supports the result obtained above tha, the interior source is substantially 
stronger (- factor of 4) than the Titanogenic source. Indeed calculations with an 
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interior source strength not substantially greater than the Titan source strength 
failed to generate an enhanced brightness above Saturn. 

4. CONCLUSIONS 


A model for the distribution of hydrogen in the Saturn system was applied to 
the best published observational data of Shemansky and Hall (1992). While the 
model for the source of hydrogen could be made more complex, it seems that the 
H Lyman-a brightness data collected by Voyager 1 and 2 for the Saturn 
encounter, as presented in Shemansky and Hall (1992), does not permit a model 
data comparison which is particularly sensitive to either the detailed spatial 
location or energetics of the hydrogen interior source. This means that a more 
sophisticated model would not be likely more illuminating. There is, however, 
current work in progress to further reduce the Voyager observations (Herbert 
1997), and it is anticipated that this will permit a more incisive model data 
comparison. Despite these limitations, it was possible to obtain very reasonable 
agreement between model and data which permitted a number of conclusions. 
Firstly, there is no question that a Titanogenic source alone is unable to account 
for the Lyman-a brightness data collected by Voyager. While a Titan source of 
3.3xl0 27 to 4.8xl0 27 hydrogen atoms s' 1 , consistent with previous estimates, was 
found to be necessary, an additional and larger by a factor of ~4 interior source of 
1.4 to 1.9xl0 28 hydrogen atoms-s' 1 was also required. This latter source implies a 
source rate of H 2 0 large compared to estimates based on sputtering of the surfaces 
in the Saturnian system, but possibly consistent with recent measurements of 
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OH. The spatial location of the tnterlor source was no, well characterized except 
that ,t must start to become large in the neighborhood of about 9 R s . The 
characteristic energy of the source hydrogens was also not well characterized 
except that the typical energy of the hydrogen from the interior source should no, 
be much below about 0.5 eV. Possibly due to the large etror bars in the data near 
Saturn, we were unable to detect any important source of hydrogen from Satum. 
Finally, it also appears that an estimate of the time averaged lifetime of hydrogen 
in the solar wind plasma at Saturn's heliocentric position is 5-7x 1 0 7 s. 
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FIGURE LEGENDS 


F. 8 ure 1. Plasma Loss Rate For H. The loss rale model of H. used in this 
study, from the plasma processes of charge exchange and electron impact 
ionization in the Saturnian magnetospheric plasma is contoured as a function of 
radial distance along Saturn's equatorial plane and vertical distance above it. The 
model is based on the Saturnian magnetospheric model of Richardson (1995 1996) 
and Richardson and Sittler (1990). The contours are in units of ,0" The 

region between - 1 2 and 20 R s with contour values of 2x 10-7 s-l actually contains 

an oscillatory lifetime structure with 2x10-7 c-l u* . , 

^ XI ° s be,n g roughly the size of the 

maxima. 

Figure 2. Comparison of Model and Data Brightness (Titanogenic 

Hydrogen only,. The Lyma„-„ brightness distribution coiresponding to the 

oyager 2 UVS scans conducted during the period 1981 DOY 180 to 186 (see 

Shemansky and Hal, ,992, is shown. The crosses are the data and solid line is the 

nrodel brightness distribution for a Titan source only with an H atom source rate 
of 4.8x10 S. 


Ftgnre 3. Comparison of Model and Data Brightness For the Period ,98, 

HOY ,80 to ,86. The Lyma„-« brightness for the Voyager 2 scans for the period 
HOY 180 to 186 in 198, are the crosses and the mode, is die solid cuive which 
consists of a Titan source of 4.8x10- atoms s ' and an interior source of I 9x10* 
atoms s a, 8.7 R, with a 40000 K Maxwellian velocity distribution. 
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Figure 4. Comparison of Model and Data Brightness For the Period 1980 
DOY 239 to 256. The Lyman-a brightness for Voyager 1 scans for the period 
DOY 239 to 256 in 1980 are the crosses while the solid line is the model 
brightness for a Titan source of 3.3xl0 27 atoms s' 1 and an interior source of 
1.4xl0 28 atoms s' 1 with a 40000 K Maxwellian velocity distribution. 

Figure 5. Contour plot of Lyman-a Brightness. The relative brightness of the 
same two source rates used in Figure 4 is contoured for the Voyager 1 observing 
geometry during the period 1980 DOY 324-343. The black dot in the center is 
Saturn and the dotted line is Titan's orbit. The brightness distribution has a 
general overall resemblance to the image in Shemansky and Hall (1992) and, in 
particular, is similarly somewhat brighter on the upper half. The Sun is to the 
right as shown by the arrow. 

Figure 6. Time Evolution of an H Atom Lost From Titan. A typical orbit for 
an H atom from Titan is displayed at various times. The orbit which evolves 
under the action of radiation pressure and the J 2 component of Saturn's gravity 
undergoes an evolution characterized by a decrease in eccentricity and a rotation of 
the major axis of the orbit. The particle eventually crashes into Saturn after 
~10 8 s. The star is towards the Sun and the numerals 2, 3, and 4 correspond to 
times of 3x1 0 7 , 5xl0 7 , and 7xl0 7 seconds respectively. 
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TABLE 1 


Physical Parameters for Titan and Saturn 


— 

Titan Exospheric Temperature (K) 

186 


Titan Radius (km) 

2575 


Titan Exobase Radius (km) 

4175 


Titan Lagrange Radius (km) 

52393 


Saturn Radius (equatorial) (km) 

60330 

— . 

Saturn Spin Rate (rad-sec' 1 ) 

1.662x10"* 


Saturn Obliquity (deg) 

26.73 


h 

1.6298xl0- 2 


Radiation .Acceleration (cm-sec 2 ) 

6.29x1 O' 3 
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DISTANCE FROM SATURN"S CENTER (Rs) 

Figure 1 
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